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Absorption spectra, exciton energy levels and wave functions for solid Ne and Ar have been calculated from 
first principles using many-body techniques. Electronic band structures of Ne and Ar were calculated using 
the GW approximation. Exciton states were calculated by diagonalizing an exciton Hamiltonian derived from 
the particle-hole Green function, whose equation of motion is the Bethe-Salpeter equation. Singlet and triplet 
exciton series up to n = 5 for Ne and — 3 for Ar were obtained. Binding energies and longitudinal-transverse 
splittings of 71 = 1 excitons are in excellent agreement with experiment. Plots of correlated electron-hole wave 
functions show that the electron-hole complex is delocalised over roughly 7 a.u. in solid Ar. 
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I. INTRODUCTION 

Electronic and optical properties of rare gas solids (RGS) 
Ne and Ar have been studied experimentaUyii2i2ii^i^i2i2i2ii£ and 
theoreticall}iii*iSii^'iiii^iiLi2ii& and have been the subject of 
several review articlesiSiSSi^. Optical absorption spectra are 
characterized by sharp exciton peaks at energies up to sev- 
eral eV below the fundamental band gap. Excitons in RGS 
consist of a hole in a p-type valence band and an electron 
in an s-type conduction band. The momentum of the hole 
can be either j = 3/2 (spin triplet) or j = 1/2 (spin sin- 
glet). Spin-orbit coupling mixes these states and so pairs of 
transverse exciton lines, labelled by principal exciton quan- 
tum numbers n and n', are observed in optical absorption ex- 
periments. When spin-orbit coupling is weak, as in Ne, ener- 
gies of these excitons ought to be close to the spin singlet and 
triplet energies. Creation of longitudinal excitons in electron 
energy loss experiments produces a longitudinal polarisation 
field which couples to the associated long-range macroscopic 
electric field and leads to energy splitting of longitudinal (L) 
and transverse (T) excitons of the same principal quantum 
number 

An integral equation approach has been applied to excitons 
in insulators for over 40 yearsa^. This incorporates screened 
electron-hole attraction and exciton exchange terms in the 
Hamiltonian. Early applications of this approach to RGSia*i^ 
used an approximation in which the exciton wavefunction 
was restricted to a single unit cell (one site approximation). 
Following this, Resca and coworkers developed an approach 
which took into account the delocalised nature of the exci- 
ton wavefunction by solving an effective mass equationiiii^ii^. 
In this approach the electron-hole attraction term was un- 
screened when both electron and hole were on a single site 
and screened by a macroscopic dielectric constant factor when 
they were on different sites. Modification of the electron- 
hole attraction term as a function of electron-hole separation 
in this way leads naturally to a quantum defect correction. 
En = Eg — Bex/{n + (5n)^, to the Wannier formula for ex- 
citon energies. En = Eg — Bex/n^- Eg is the fundamen- 
tal band gap and B^x is the exciton binding energy. Resca 
and Restai^ showed that the former expression could predict 
exciton energies rather well with a weak dependence of the 



quantum defect, on n. However, subsequent direct mea- 
surements of the fundamental band gaps in RGS~- showed 
that the fundamental gap derived from a fit to the Wannier 
formula (excluding the n' = 1 exciton) yielded a value for 
Eg in good agreement with the experimental values while the 
quantum defect model yielded quite different values. Bern- 
storff et alrf^ concluded that RGS excitons 'do not possess 
atomic parentage'. Exciton wavefunctions for RGS are there- 
fore well known to have strong but incomplete localization 
of the electron and hole on the same site. This intermedi- 
ate binding character of excitons in RGS means that they are 
far from the well-separated electron-hole pair limit, which 
applies to semiconductors and is well described by the ef- 
fective mass approximation theorya^. A recent calculationii, 
which used a screened Coulomb electron-hole interaction and 
a Slater-Koster parametrization of the band structure, found 
good agreement with experimental exciton energies and delo- 
calization of the electron-hole wavefunction over three nearest 
neighbor distances. 

The integral equation approach was applied to diamond2^i2i 
and silicor^S by Hanke and Sham in the 1970's using a tight- 
binding parameterization of the band structure. Strinati2S 
showed how the equation of motion for the particle-hole 
Green function, the Bethe-Salpeter equation^ for core-hole 
excitons could be reduced to an effective eigenvalue prob- 
lem. Recent ab initio calculations of valence excitons in crys- 
talline solids^ii^*iS have been based on this effective eigen- 
value problem and the Bethe-Salpeter formalism has been 
reviewed^. In this paper we present ab initio calculations 
of excitonic absorption spectra and correlated electron-hole 
wave functions for excitons in Ne and Ar. Quasiparticle 
band energies are calculated using the GW approximation^ 
and the Bethe-Salpeter formalism^, which includes statically 
screened electron-hole attraction and exciton exchange terms, 
is used to calculate the optical spectrum. This is generated us- 
ing matrix elements between the ground state and correlated 
electron-hole excited states and the longitudinal-transverse 
(LT) splitting of excitons is investigated. 

The remainder of the paper is arranged as follows: In Sec- 
tion|n]the theoretical formalism is presented, in Sectionlllllre- 
sults of calculations of optical spectra and correlated electron- 
hole wave functions in solid Ne and Ar are given and finally 
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conclusions are given in Section HVl 



II. THEORETICAL BACKGROUND 

A. Quasiparticle energies 

The starting point in our approach is to generate the quasi- 
particle energies and wave functions of the system. Quasi- 
particle energies are obtained by solving the quasiparticle 



equatioi 
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i/(r)^«^(r)+ J S(r,r',i?)V'r(r')dr' 



(1) 

using perturbation theory. The unperturbed Hamiltonian is a 
Kohn-Sham Hamiltonian and the self-energy operator, E, is 
obtained using the GW approximation (GWA). Quasiparti- 
cle wave functions, ■0j^^(r), are approximated by eigenfunc- 
tions of the unperturbed Kohn-Sham Hamiltonian. Density 
functional theory (DFT) within the Perdew-Wang generalized 
gradient approximation^ (PWGGA) is used to obtain DFT 
eigenvectors and eigenvalues. Details of GW calculations as 
well as quasiparticle band structures along symmetry lines for 
Ar and Ne are given in Ref. IH. The CRYSTAL coda^ 
was used to generate single-particle wave functions for Ne and 
Ar in an all-electron Gaussian orbital basis and the Coulomb 
potential was expanded in plane waves. GW and BSE cal- 
culations were carried out using the EXCITON^ code. The 
spin-orbit interaction was omitted from the calculations and 
experimental lattice constants were usedi^. 



B. Electron-hole excitations and optical spectra 

Correlated electron-hole states, \N, S), can be represented 
in a basis of single-particle quasi-electron (conduction, c) and 
quasi-hole (valence, v) states. 



\N,S) 



t6t|7V,o)=.^ALc|k^;c}, (2) 



where and fo^ create a quasi-hole and quasi-electron, re- 
spectively, in the many -body ground state |A^, 0). Coupling 
coefficients, A^^^, and excitation energies, ils, are obtained 
by solving the BSE in the form™, 



lz.'v' c' 



S 

live 1 



(3) 



where 5 denotes the electron-hole interaction and energies 
i?kc and iJkt) are quasiparticle energies obtained within the 
GWA. The interaction kernel, S, is given as a sum of two 
termg-4: the screened electron-hole attraction, also called the 
direct interaction, S'', and the exchange interaction, S"^, which 
results from bare Coulomb repulsion. Neglecting any dynami- 
cal screening, the matrix element of the direct term in a plane- 



wave basis is given by, 

^ ' ^~ ^ ^^Jq+G||q + G'| 

x(w'k'|e'(i+^>'-|«k)(ck|e-*(i+<=^')"-|c'k')(5q^k'-k. (4) 

£qq, (q, 0) is the symmetrized, inverted, static RPA di- 
electric matrijiffi, G and G' span the reciprocal lattice and 
denotes the crystal volume. When computing the matrix el- 
ements in Eq. Q special care has to be taken for the case 
q — > 0. If G = G' = the interaction diverges as Xjq^ . This 
contribution is separated from the left side of Eq. (|3} and in- 
tegrated over a small sphere of volume V = Vbz l^k where 
Vbz is the volume of the Brillouin zone and A^^ is the num- 
ber of k points, as suggested by Arnaud and Alouani-^. This 
divergence contributes notably only when v = v' and c = c'. 
In addition a divergence of type 1/q occurs when one of the 
G vectors is zero. These terms are neglected, because their 
contribution either averages to zero or vanishes in the limit of 
a large number of k points^. 

The exchange term of the interaction for singlet states has 
the form^ 



(wck|S^|i/c'k' 
47re^ 



2 X 



^ij(ck|e^«--|«k)(«'k'|e-«-"-|c'k'). 



G/0 



(5) 



The G = term omitted from Eq. (|5} is responsible for LT 
splitting of singlet excitons^. It has the form, 



2 X 



Aire' {ck\p\vk) ^^j, (^^^k^lplc^kQ 

^ Eke — £^ku ^'k'c' — Ei^'yi 



(6) 



where Q = q/|q| is a unit vector parallel to the macroscopic 
field and p is the momentum operator It is obtained by re- 
placing G by q in Eq. |6j taking the limit q ^ and using 
the commutation relation [H, r] = ip to avoid calculating 
matrix elements of the undefined dipole operator. It has been 
omitted from recent ab initio calculations of the transverse 
(optical) spectrum in semiconductors^ifS. The LT splitting is 
generally small (< 0.1 eV) in that case. However in Ne and 
Ar it is large enough to warrant its inclusion, when ab initio 
results are to be compared to experiment. The exchange term 
in the interaction vanishes for triplet exciton states. 

Finally, Eq. (|3} is solved and the macroscopic dielectric 
function is given hy^ 



SAiiu) = 1 + lim 7-^ 



E 



E.,c.k(^k|e-'q-|ck)Ag„J 

ns-uj- iQ+ 



(7) 

When the term in Eq. (|6} is omitted the LT splitting induced 
by exchange terms in Eq. |5]with G is around 0.02 eV, but 
when it is included these quasi-degenerate modes split into L 
and T modes with optical transition moments either parallel 
(L) or perpendicular (T) to the unit vector Q. The splitting is 
large (> 0.1 eV) only for n' ~ 1 singlet excitons. 
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In Eq. Q optical transitions are given as a coherent sum 
of the transition matrix elements of the contributing electron- 
hole pair configurations, including the coupling coefficients, 
^kuc- Without the electron-hole interaction, excitations are 
given by vertical transitions between independent electron and 
hole states. In that limit Eq. reduces to the well known 
RPA dielectric function. 
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.RPA 
'M 



(w) = 1 + 2 lim 



87re^ 



(wk|e-*i-''|ck)|2 



(8) 

In calculating £m from Eq. 0, three valence bands, one 
conduction band and 2048 Monkhorst-Pack'*^ special k points 
in the full Brillouin zone were used. An artificial broad- 
ening of 0.05 eV was introduced in spectral lines. Several 
groups^ii^i^i have used low symmetry, shifted k points to 
achieve well converged excitonic spectra. Excitons in RGS 
belong to an intermediate regime, where they are localised in 
real space and delocalised in reciprocal space (see Section lini i. 
Convergence of excitonic spectra can therefore be achieved 
using special points only. 65 G vectors were used to calculate 
the direct part of the interaction, Eq. (0}, and up to 500 were 
used for the exchange interaction, Eq. (|5}- An all-electron 
basis containing 5 s, Ap and 2d Cartesian Gaussian orbitals on 
the atomic nuclear site was used for Ne and a similar basis 
with 7s, 6p and Ad orbitals on the atomic nuclear site plus 2s, 
2p and Id orbital on the octahedral interstitial site was used 
for Ar. The basis used for Ne is smaller than that used pre- 
viously for a GW calculation on Ne while the Ar basis is the 
same as used previously (basis set 2 in Ref. |35||). 



III. RESULTS 
A. Optical Spectra 

1. Neon 

The BSE eigenvalue problem was solved for Ne using an 
exciton exchange term which either included or omitted the 
LT splitting term (Eq. (|6}). Exciton energies and binding en- 
ergies, Eb — Eg — En, are given in TableUand are compared 
to earlier calculations. In the work by Andreoni et al^, n = 1 
and n' = 1 excitons were calculated using matrix elements 
and a band structure from an augmented plane wave (APW) 
calculation and excitons with higher principal quantum num- 
bers were calculated using the effective mass approximation 
(EM A). In the work by Martinelli and Pastori Parravicinii^ 
the electron-hole attraction term in the integral equation was 
treated using a model screened Coulomb potential. Experi- 
mental exciton energies in Table U are from optical transmis- 
sion data by Saile and Koch^. Both L and T excitons were 
observed in thin film optical absorption data and the L exci- 
ton energy of 17.75 eV from the BSE calculation is in good 
agreement with a value of 17.74 eV from electron energy loss 
datiO^. 

The imaginary part of the macroscopic dielectric function, 
62, derived from Kramers-Kronig transformation of reflec- 
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FIG. 1 : Imaginary part of the macroscopic dielectric function for Ne 
calculated using BSE wave functions and from experimental data by 
Skibowski— . 



TABLE I: Exciton energy levels, En, band gaps. Eg, binding en- 
ergies, Eb and LT splittings in solid Ne in eV. Theoretical results 
from an augmented plane wave (APW) calculation for the n — 1 
and n' = 1 states and an effective mass approximation (EMA) cal- 
culation for n and n' > 1, a model screened electron-hole poten- 
tial calculation and a BSE calculation are compared to experimental 
peak positions in optical transmission data. The APW and EMA 
calculations include spin-orbit coupling. Fundamental band gaps 
and longitudinal-transverse splittings, ALT, are given in the last two 
rows. 



n 


En" 


En'' 


En 


Eb' 


En'' 


Eb' 


1 


17.51 




17.25 


4.44 


17.36 


4.22 


r 


17.86 


17.37 


17.45 


4.24 


17.50 


4.08 


2 






19.90 


1.79 


20.25 


1.33 


2' 


19.98 


20.64 


19.95 


1.74 


20.36 


1.22 


3 






20.55 


1.14 


20.94 


0.64 


3' 


20.93 


21.19 


20.55 


1.14 


21.02 


0.56 


4 






20.95 


0.74 


21.19 


0.39 


4' 


21.25 


21.40 


20.95 


0.74 






5 






21.15 


0.54 


21.32 


0.26 


5' 


21.40 


21.50 


21.15 


0.54 






E, 


21.67 


21.69 


21.69 




21.58 




ALT" 


0.23 




0.30 




0.25 





"APW/EMA calculation Ref.Q 
*SK calculation Ref.jil 
'BSE calculation This work. 
''ExperimentRef.il 
' j = 1/2 exciton. 



tion measurements by SkibowskiiS and a BSE calculation (Eq. 
(0) for singlet states without the LT splitting in the exchange 
term are shown in Fig. ^ The experimental spectrum shows 
excitonic absorption at 17.49, 20.24, 20.87 and 21.30 eVi^. 
These are the n' ~ 1, 2, 3, 4 (singlet) exciton energies for Ne 
and they are nearly all coincident with data for Ne by Saile 
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and Kochi (Table U, where data for triplet exciton energies 
are available also. Singlet and triplet BSE calculations show 
the first two exciton absorption features at 17.45(17.25) and 
19.95(19.90) eV. While there is good agreement between re- 
sults of the BSE calculation and experiment for n' and n = 1 
states, exciton binding energies for n' and n > 1 are overes- 
timated. This may be a result of incomplete screening of the 
electron-hole interaction caused by the limited basis set used 
in this calculation. Differences in n and n' binding energies 
decrease with increasing quantum number Singlet/triplet en- 
ergy splittings in BSE calculations for the first two states are 
0.20 and 0.05 eV and these compare to 0.14 and 0.11 eV in 
experiment (Table HJi. For higher states the BSE predicts es- 
sentially no splitting while a small splitting is still found in 
experiment. 

When the LT splitting term (Eq. |6} is included in the singlet 
state exchange, the n' = 1 exciton shifts to 17.75 eV, which 
corresponds to an LT splitting of 0.30 eV. The experimental 
LT splitting energy in the n = 1' exciton is 0.25 eV and the 
calculation by Andreoni et al^ gave a value of 0.23 eV. The L 
exciton energy coincides with the experimental value of 17.75 
eV determined by optical absorption in a thin filrrA and a value 
of 17.74 eV obtained from electron energy loss experimentsiS. 
The LT splitting for the n' = 2 exciton from a BSE calculation 
is much smaller and has a value of 0.03 eV. Energies of T 
excitons remain essentially the same when the exchange term 
in Eq. is included. 



2. Argon 

Singlet and triplet exciton energies and binding energies for 
Ar derived from BSE calculations and experiment are given 
in Table |ll] Experimental exciton binding energies in Ar are 
reproduced well in parametrized calculations by Andreoni et 
alii^ and by Baroni et alM-. Binding energies from BSE calcu- 
lations are also in good agreement with experiment, although 
BSE exciton energy levels lie below experimental values be- 
cause the band gap is underestimated by 0.5 eV in the GW 
calculation used to generate the quasiparticle energies in Eq. 
|3l Exciton binding constants for j = 3/2 and j = 1/2 states 
have been estimated from experimental^. By fitting n and 
n' > 1 exciton peak positions to a Wannier plot (i.e. with no 
quantum defect) we find B^x to range between 2.06 and 2.23 
eV for the j = 1/2 exciton and between 2.17 and 2.30 eV for 
the j ~ 3/2 exciton. These values are in good agreement with 
collected data by Bernstorff et al^. When the n and n' = 2 
levels from the BSE calculation are fitted we find a binding 
constant of 2.56 eV and when we use the n and n' = 3 lev- 
els we find a binding constant of 2.16 eV. The experimental 
LT splitting energy in the n = V exciton is 0.15 eV. The 
LT splittings derived from the BSE calculation (0.36 eV) is 
significantly larger than that derived from the calculation by 
Andreoni et al^ (0. 19 eV) and the experimental value. 

The imaginary part of the macroscopic dielectric function 
from experiment and a BSE calculation for the j = 1/2 exci- 
ton are shown in Fig. |2l The experimental optical spectrum 
contains both the n = 1 and n' = 1 peaks because spin- 



TABLE II: Exciton energy levels, En, band gaps. Eg, binding en- 
ergies, Eb, LT splittings and binding constants, Bex, in solid Ar in 
eV. Theoretical results from an augmented plane wave (APW) cal- 
culation for the n = 1 and n' = 1 states and an effective mass ap- 
proximation (EMA) calculation for n and n' > 1, a model screened 
electron-hole potential calculation and a BSE calculation are com- 
pared to experimental peak positions in optical transmission data. 
The APW and EMA calculations include spin-orbit coupling. Funda- 
mental band gaps, longitudinal-transverse splittings, ALT, and bind- 
ing constants from are given in the last three rows. 



n ±!jb J^b 


Jin 


Eb 


TP d 


TP d 

Eb 


1 2.12 


11.60 


2.09 


12.10 


2.06 


1' 1.84 1.86 


11.75 


1.94 


12.35 


1.90 


2 


13.05 


0.64 


13.58 


0.58 


2' 0.47 


13.05 


0.64 


13.75 


0.50 


3 


13.45 


0.24 


13.90 


0.26 


3' 0.19 


13.45 


0.24 


14.03 


0.22 


E, 


13.69 




14.25 




ALT 0.19 


0.36 




0.15 




Bex 


2.56 




2.06 




"Refill 










*Ref.Il2 










'^BSE calculation This work. 










''Experiment Ref . |2 










' j = 1/2 exciton. 










TABLE III: Fitted exciton binding constants, B^x, and band gaps. 


Eg, in solid Ar in eV. n and 


n' > 1 levels have been used for raw 


experimental data. The n = 


2 and ; 


n' = 2 enerj 


;ies and the GW 


value for Eg have been used for BSE data and result in 


the same 


binding energies for j = 1/2 and j = 


: 3/2 series. 






j Bex" Eg" 


R * 


E/ 


Bex 


E/ 


111 2.06 14.25 


2.23 


14.31 


2.56 


13.69 


3/2 2.30 14.16 


2.17 


14.12 


2.56 


13.69 



"Ref. 13 
"Ref. lis 
"This work. 



orbit coupling mixes singlet and triplet states thereby creating 
two optically active transitions. The GW band structure used 
here did not include spin-orbit coupling and the BSE spectrum 
shown includes only optically active singlet excitons. 

B. Electron-hole wave functions 

The ab initio approach used here allows the wave function 
of an electron-hole excitation to be examined in detail. Ex- 
cited state wave functions, Xn{^e, ^h), in a coordinate repre- 
sentation are, 

occ empty 

Xnire,r^) = EE E AZ^cMeWkM. (9) 

k u c 

The coordinates and refer to the position of the electron 
and hole, respectively. The wave function is a six-dimensional 
scalar function. In a crystalline system, it is invariant to lattice 
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FIG. 2: Imaginary part of the macroscopic dielectric function for Ar 
calculated using BSE wave functions and from experimental data by 
Saile— . The fundamental gap predicted by the GWA calculation for 
Ar is less than the experimental gap by 0.5 eV. The scale for the BSE 
spectrum has been shifted by 0.5 eV in order to align fundamental 
gaps in theory and experiment and facilitate comparison of experi- 
mental and theoretical spectra. 
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FIG. 4: Same as Fig.|3|along the line AB: electron distribution (left 
panel) and hole distribution (right panel). Blue circles present atom 
positions. 



onto nearest neighbor sites. A more quantitative display of 
the electron-hole correlation function is obtained by plotting 
the modulus squared of the wave function along a line con- 
taining three Ar atoms (Fig.0). These results confirm ear- 
lier conclusions^iii*^ that the first exciton is delocalized over 
nearest neighbor atoms and give for the first time detailed in- 
formation on the structure of correlated electron-hole func- 
tions in Ar in real space. 




FIG. 3: Real space probability density (|x(rhi I'e)!'^) for an electron 
(re) with respect to a fixed hole (r^) (left panel) and the distribution 
of a hole with respect to a fixed electron (right panel) in the (100) 
plane of solid Ar for the n' — 1 state. The hole (electron) is fixed 
at the central atom (blue circle). Probability densities along the line 
AB are presented in Fig. |4] Green empty circles correspond to atom 
positions. 



translations simultaneously applied to Vh and Ve- By fixing ei- 
ther the electron or hole in space and plotting \xn{^e', r/i)P, 
details of the spatial correlation function for the electron-hole 
distribution can be explored. Fig. |3] illustrates the probabil- 
ity density for an electron when a hole is fixed at an atomic 
nuclear site as well as the probability density for an elec- 
tron with a hole fixed at the atomic position. The plot with 
the fixed hole shows an electron distribution in predominantly 
atomic s states with a p lobe distribution over atomic sites. 
The plot with the electron fixed shows the hole to be pre- 
dominantly on the same atomic site with some delocalisation 



IV. CONCLUSIONS 

An exciton Hamiltonian has been diagonalized to obtain 
singlet and triplet exciton series for Ne up to ?i = 5 and Ar 
up to n = 3. Exciton binding energies for n' = 1 singlet 
and n = 1 triplet excitons for both Ne and Ar are in excellent 
agreement with experiment. The longitudinal-transverse split- 
ting of the ji' = 1 singlet excitons in Ne is in good agreement 
with experiment; while the value obtained from a BSE cal- 
culation for Ar exceeds the experimental value (0.36 eV c.f. 
0.15 eV), one might expect Ar to have a larger splitting than 
Ne, as predicted by the BSE calculation (0.36 eV c.f. 0.30 
eV), because of a larger polarisability density in Ar. 

The band gap for Ne derived from a GW calculation is 
21.69 eV and agrees very well with the experimental value 
of 21.58 eV. The band gap for Ar from a GW calculation is 
13.69 eV and underestimates the experimental value of 14.25 
eV by 4%. The excellent agreement between theory and ex- 
periment for the band gap of Ne is fortuitous and is a result of 
using a relatively small basis set with functions located only 
on atomic sites. When s and p Gaussian orbitals were added 
to the basis at octahedral interstitial sites of the fee lattice, 
the GW band gap was reduced to 20.04 eVi^ an underes- 
timate of the experimental value by 7%. Underestimation of 
experimental band gaps to this extent is typical of perturbative 
GW calculations which start with self-consistent Kohn-Sham 
Hamiltonians, as this work does. 

Overestimation of exciton binding energies with n > 2 for 
Ne, but not for Ar, for which a superior basis set was used, 
suggests that the limited basis set used for Ne does not ade- 
quately account for screening of the electron-hole attraction 
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for more extended excitons. When a BSE calculation was 
performed for Ne using a basis which contained additional 
interstitial site basis functions, the optical spectrum contained 
spurious features above the first exciton absorption, although 
the GW calculation with this basis set-'^ did not suffer from 
similar problems. Binding energies of n = 1 and n' ~ 1 exci- 
tons are expected to be much less affected by underscreening 
than those with larger principal quantum numbers since the 
electron and hole are in close proximity in those states and 
Coulomb potential is poorly screened at that range. 

Plots of correlated electron-hole wave functions for the 
n = 1 exciton in Ar show that the extent of delocalisation 



of the electron-hole pair is quite limited, confirming that these 
excitons belong to an intermediate regime between Frenkel 
and Wannier types. 
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